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Abstract 



Using perturbative renormalization group we study the influence of ran- 
dom velocity field on the critical behaviour of directed bond percolation 
process near its second-order phase transition between absorbing and active 
phase. We consider Kraichnan model with finite correlation time for mod- 
elling advecting velocity field. Using functional integral representation we 
are able to apply field-theoretic renormalization group to determine possible 
universality classes. The model is analyzed near its critical dimension by 
means of three-parameter expansion in e,5,r], where e is the deviation from 
the Kolmogorov scaling, 5 the deviation from the critical space dimension 

^jI dc and rj is the deviation from the parabolic dispersion law for the velocity 

• correlator. Fixed points with corresponding regions of stability are evaluated 

^ j to the leading order in the perturbation scheme. 

T! 1 Introduction 

> 

Directed bond percolation (DP) problem is one of the most famous model in sta- 
tistical physics exhibiting non-equilibrium second-order phase transition [T]. In 
C^ various formulations it can serve for explaining hadron interactions at very high 

energies (Reggeon field theory) [51 [3] , various models of disease spreading [3 H] 
or as in original formulation f5j wetting of porous material or exploring path in 
labyrinth. The upper critical dimension for this problem was estimated to be 
d = 4 in contrast to the value rf = 6 for the isotropic dynamical case [2 El E] j 
which we do not consider. It was conjectured [H E] that critical regime of any 
one-component system with short-range interactions with continuos transition to 
absorbing state belong always to the DP class. Despite the general validity of this 
statement the percolation model was not experimentalluy observed until recently 
|10| . In this work transition between two topologically different states in elec- 
trohydrodynamic convection was observed. Diffuculties with direct experimental 
observation of percolation results from the fact that real material usually contains 
inhomogenities, various defects, anisotropics and so on. They can destroy or make 
observation of phase transition completely hopeless. Therefore it can be of im- 
portance to study possible deviations from the simple model of percolation and 
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quantify their effects on on percolation process. A lot of effort was put into the 
investigation of various effects, e.g. in papers |lll [T^ [TB] long-range interactions 
by the means of Levy-flight jumps were studied both in time and space variables, 
introduction of immunization was examined in [H [TSl [IJ , effect of surfaces was 
studied in [15J etc. One can also easily imagine that spreading of disease can be 
rapidly enhanced by some external atmospheric current or by flying insects. In 
both cases because of additional drift can be modelled as random velocity field j T7] 
with prescribed statistical properties. In this paper the influence of advective fleld 
described by rapid-change Kraichnan model was studied, which is characterised by 
white-in-time nature of velocity correlator. Generalizing this approach it is pos- 
sible to study e.g. effect of compressibility [18, or using stochastic Navier-Stokes 
equations effects of "real" turbulent field [12]. In this work we try to investigate 
the infiuence of finite correlated velocity field (for introduction see [5D]) and de- 
termine how it can change the critical behaviour of percolation process. By the 
means of renormalization group approach we determine possible fixed points with 
corresponding regions of stability. We show that the model exhibit 10 possible 
large-scale regimes. 

This paper is organized as follows. In section 2 we give the detailed description 
of the model. In section 3 we perform dimensional analysis (power counting) of 
the model and prove its multiplicative renormalizibility. In section 4 we briefly 
describe main ingredients of the diagrammatic technique and present flrst order 
(one-loop) calculation of the renormalization constants, which are given in the ex- 
plicit form. In section 6 we analyse asymptotic behavior of the model according 
to its fixed point structure. We give a comprehensive account of corresponding 
critical exponents and present range of stability in the (e, S, rj) space. Section 6 is 
devoted for conclusions and future plans. 

2 Field theoretic formulation 

The most rigorous approach to the percolation problem is based on the interpre- 
tation of it in the means of reaction-diffusion process. According to the standard 
approach[T] master equation for such problem can be rewritten employing Doi 
formalism|21| into the form of time-dependent Schrodinger equation with non- 
hermitean hamiltonian. After performing continuum limit the effective action can 
be derived, which is amenable to the usual field-theoretical methods. However it 
can be shown[T] that the phenomenological approach based on the use of Langevin 
equation with suitable chosen noise leads to the same prediction of universal quan- 
tities as the aforementioned master equation approach. Let us therefore briefly 
describe the main points of the latter. The stochastic non-linear differential equa- 
tion for the coarse-grained density of infected individuals(agents) ip{t,x.) can be 
written in the following form jT] [^ 

dMt, x) = DoiW' ~ ToMt, x) - ^i^^t, x) + at, x), (1) 

where V^ is Laplace operator. Do is a diffusion constant and Ao is a positive 
coupling constant. The parameter tq measures the deviation from the threshold 



value for the infection probability. One can assume tq '^ Pc — P, where p^ is 
a critical probability for observing percolation (analogous to the deviation from 
critical temperature for equilibrium models). It is important to note, that the 
model has unique absorbing state with 'i/;(i,x) — (no sick agents) from which it 
cannot escape. The Gaussian short-ranged noise Ci^) with zero mean accounts for 
the density fluctuations and it has to respect the absorbing state condition. This 
can be achieved by the following choice of its correlatorjl] [TTJ [T7] 

(C(i, x)C(i', x')) = i?oAoV(i, x)5(i - t')<5(x - x'). (2) 

In this work we want to apply renormalization group technique for the study of 
large scale behavior of the model (jTl. It is therefore useful to recast it into the 
path integral formulation. The Langevin equation (nj) has the standard form of 
stochastic dynamic problem|23]. By introducing Martin-Siggia-Rose|24| response 
field V'^(i,x) and integrating out the Gaussian noise it is possible to obtain action 
functionaHS) i24j for the pure directed percolation problem or Reggeon field theory 
respectively 

5lo(V'^V') = V^(-a,+i?oV^-i?oro)^ + ^[(V'^)V-V'V^']. (3) 

For convenience the required integrations over the space-time variables are not 
explicitly indicated in the action (Isl, e.g. the second term means the following 
expression 

V^ty^V'^ /rfi /rfx V^(t,x)V2V'(t,x). (4) 

The model (pi) satisfied the so-called rapidity reversal symmetry (in the language 
of Reggeon field theory) 

V'(i,x)^-Vt(_t^x), V^(i,x)^-^(-f,x), (5) 

which should be respected by the renormalization transformation. 
In this work we would like to study the influence of advective flcld on the spreading 
of agents ip. The agents can be considered as passive scalar quantity j^HI that is 
advected by the velocity fleld with no back influence on the velocity fleld itself 
with nontrivial interactions given by the cubic terms in ([3| . The inclusion of the 
velocity fleld v(i,x) corresponds to the replacement 

dt-^Wt^dt + (v.V), (6) 

in equations M and (pi), where Vt is convective (Lagrangian) derivative. The 
most realistic description is based on the use of stochastic Navier-Stokes equations 
|23| I26| . However, in this paper we shall study a simplifled model in which we 
prescribe statistical properties of the velocity fleld. Let assume that v{x) is a 
random Gaussian variable with zero mean and the correlator given as 
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Here Pij(k) = Sij — kikj/k"^ is the transverse projection operator, k — |k| is the 
wave number and the kernel function D^ is assumed to have the following form 

^, ,, gi,Dlk^~^'-^^-^ ,,, 

^''^"'^) = C.2+,2^2(fc.-,). - (8) 

Here gio is the coupling constant (small parameter of the ordinary perturbation 
theory) and the exponents e,^ and 77 play the role of small expansion parameters. 
They could be regarded as an analog of the expansion parameter e = 4 — d in 
the usual sense of dimensional regularization. However, in this paper e should be 
understood as deviation of exponent of the power law from that of the Kolmogorov 
scaling \^ , whereas 6 is defined as the deviation from the space dimension two via 
relation d = A — 25. The exponent 77 is related to the reciprocal of the correlation 
time at the wave number k. The parameter uq may be used for labelling of the 
fixed points and has the meaning of the ratio of velocity correlation time and the 
scalar turnover time |28] . Although in the real calculations e is treated as a small 
parameter, the real problem corresponds to the value e = 4/3. 
It is interesting to note that the model for the advection field v(x) contains two 
cases of special interest: 

(a) in the limit uq — > oo^g'-^Q = gio/ug ~ const we get the 'the rapid-change 
model' _Di,(a;,k) — 7> g'^^Dok^^^^^^^'^^^ , which is characterized by the white- 
in-time nature of the velocity correlator. 

(b) limit uq — > 0, ^"o — .910/^0 = const corresponds to the case of a frozen velocity 
field Dy{u!, k) — ?► g"QDQTr6{uj)k^^~^'^ , when the velocity field is quenched (time- 
independent). 

The averaging procedure with respect to the velocity field v(x) may be performed 
with the aid of the following action functional 

S'20 = — I dtdt' dxdx' v{t,x)D-'^(t-t',x-x')v{t',x'), (9) 

where D~^ is the inverse correlator (|7| (in the sense of the Fourier transform). 
The expectation value of any relevant physical observable may be calculated using 
the complete weight functional W {ip^ , ip , v) — e"^i°+"^^°, where Si and S2 are the 
action functionals (l3| and (|9|. 

The full problem is equivalent to the field-theoretic model of the four fields 
$ = {ip^ , tp, v} with the total action functional given as the sum of functionals (p| 
and @ 

50($)-^10(*)+^20(3')- (10) 

Formulation ( 10 1 together with (l9| and (Isl means that statistical averages of ran- 
dom quantities can be presented as functional averages with the weight exp 6*0 ($) 
and the generating functionals of total G{A) and connected W^(A) Green functions 
are represented as functional (path) integral 

G{A) = expW{A)= /"x>$exp[S'o($) + A$] (11) 



with sources A = {Ay,Ay,A^f,A^} as a scalar product A^ — ^^v4$$. All 
correlation and response functions can now be calculated by perturbative means 
in a standard fashion [23J US]. The general consequence of causality is the vanishing 
of Green functions 

(l^^(tl,Xi)V'^(t2,X2)...V^(tw,Xjv)), (12) 

which should be satisfied for arbitrary N for any stochastic model[^ . Let us also 
note that the inclusion of the velocity field does not break the symmetry (|5| . By 
the direct inspection of perturbation theory it can also readily be seen that the real 
expansion parameter is rather Ag than single Aq. This fact can also be understand 
as a consequence of the symmetry (pi) . Therefore we introduce new charge by the 
relation 

920 = Aq, (13) 

which corresponds to the charge uq used in the literature [11 II 7|. 



3 Scaling analysis and UV renormalization proce- 
dure 

The theoretical analysis of the UV divergences is based on the power counting 
analysis[5S]. In contrast to the static models quantities in dynamical models are 
invariant under two independent scale transformations (with respect to time and 
space variable). Therefore the canonical dimension of quantity Q is fully deter- 
mined by two canonical dimensions, the frequency dimension dn and momentum 
dimension dg . These dimensions are found are found from the usual normalization 
condition 

d^ = -dr = l,dfe = -rf'-l,dfc-rf'=0. (14) 



and from the requirement that the action ( 10 1 is dimensionless with respect to the 
momentum and frequency dimensions separately. The total canonical dimension 
dq = dn + 2dQ is determined from the condition that the parabolic differential 
operator of the diffusion Dq'V^ and time differential operator dt scale uniformly 
under the transformation k — > /ifc, lu — > ^^lu (corresponds to the free theory with 
dt oc Dq^ ). It plays the same role as the canonical (momentum) dimension for 
static models. The canonical dimensions of the model ([To|) are given in the Table 
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dq 


d/2 


d/2 


1 


2e + 77 


5 


25 





■q 
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Table 1: Canonical dimensions of the fields and bare parameters 
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fl] We see that the model is logarithmic (the canonical dimensions of the coupling 
constants 510,320 ^^nd uq simultaneously vanish) at space dimension d = 4 (or 
equivalently (5 = 0) and for the choice e ~ 77O. In what follows we will employ 
dimensional regularization with the minimal subtraction (MS) scheme. According 
to the general theory ||25j of the renormalization group the UV divergences in the 
Green functions in this scheme manifest themselves as poles in e, S, rj or possibly 
as their linear combination. 

The total canonical dimension for arbitrary one-particle irreducible (IPI) Green 
function F = ((f> . . . $)i-ir is given by the relation [25] 

dT = d + 2- N^d^, (15) 

where N^ = {-/V^,, A'^t , N^} are the numbers of the external fields entering into the 
Green function T and the summation over all types of fields is implied. From the 
symmetry (|5| it follows that the counterms corresponding to the terms (?/'^)^'0 and 
■f/j^i/)^ can be renormalized by the same renormalization constant. We are thus led 
to conclusion that all terms that should be renormalized are already present in the 



action. The model ( 10 1 is multiplicatively renormalizable and the renormalized 



action for it can be written in the general form 

The renormalization action can be obtained by the multiplicative renormaliza- 
tion of the fields tp^ — > Z^,-tip\ tp — )• Z^fjtp and the parameters 

Dq = DZd, gioSd^^giii^^^'^Zg^, Ta=TZr, uo^ufi'^Zu, 

Ao = A^ Zx, 92oSd = 252M^ Zg^, (17) 

where Sd = Sd/{'2TT)'^ {Sd is the volume of unit sphere in d-dimension) is the 
common factor resulting from the momentum integration and factor 2 was inserted 



for the convenience. From the renormalized action (16 1 and the definition (17 1 it 
is easy to relate renormalization constants as follows 

Zi = Z.ti,Z.^,-\ = Z.tiiZ^tZy, Z2 ~ Z.^Z^i- Z]j, Z^ ~ Z^Z^i Zt-Zd , 

Z4 = ZdZxZ^p'iZjp = Z]jZ\Z^,fZ^, Z5 = Zi, (18) 



Since the nonlocal term involving v fields in ( 16 1 should not be renormalized, the 



relations Zg-^Z'^ = 1 and Z^Zu = 1 have to be satisfied. Inverting relations (18 1 
leads to the relations 

Zy = 1, Z^ =^ Z^,f ~ Z^ , Z\ = Z4Z2 Z^ , Z]j ^ Z2Z^ , 

Zu = •^l'Z^2 1 ^T = Z3Z2 , Zg^ = Z1Z2 , Zg.^ = Z^Z^ Z2 ■ (19) 

4 Calculation of the renormalization constants 

The standard perturbative approach is based on the diagrammatic expansion into 
the Feynman graphs |231[25J. The inverse matrix of the free (quadratic) part of the 
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Figure 1 : The propagators of the model 
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Figure 2: Interaction vertices of the model 

actions dsl and (pi determines the form of the bare(unrenormalized) propagators. 
In the Feynman graphs these propagators correspond to lines connecting interac- 
tion vertices. It is easy to see that the studied model contains three different types 
of propagators graphically depicted in the Fig. 1 . This has to be contrasted with 
the models [TTJ [16] where the inclusion of long-range interactions led just to the 



redefining ip^tp propagator. The propagators for the model (10 1 have the following 
form 



uj^ + ulDlik"^-^) 



<^">0 = ^.2^!,2n2a2-.^2 -P'^W 



{^i^% = {ij^i^)l^^-—-^-——- (V-V;)o = (V'V)o = (20) 
-lujk + Do[k^ + To) 

in the frequency-momentum (oj, k) representation. The vertex fact or P5] 

(5™l/($) 

is associated to each interaction vertex of Feynman graph. Here, $ could be any 
member from the full set of fields {tp^^jv}. From the action (|3| we thus derive 
three possible interaction terms depicted in the Fig. 2. First two of them are the 
usual interaction vertices for direct percolation or Reggeon field theory. The last 
one is responsible for the advection of spreading agent by the velocity field. In 
the case of incompressible fiuid it is convenient to transfer the derivation on the 
response field -0^ 

dtdxip^ 'V (vip) = / dtcbap^ di{viip) = — I dtdx^^ Vidiip = / dtdK{diilj^)vi^, 
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where we have used the fact that fields vanish in the infinity. Rewriting this 
expression in the form -ip^VjVjip we obtain immediately the vertex factor in the 
momentum space 

Vj = ikj. (22) 

From the explicit form of propagators and interaction vertices perturbation series 
in Feynamn diagrams for the (connected or irreducible) Green functions could be 
constructed. From the condition of UV finitness of IPI functions r^,t^ and r(^,t)2^ 
all renormalization constants Zi — Z4 can be calculated. One-loop approximation 
leads in MS scheme to the following result 



Zi = l 

Z3 = l 



52 

8(5' 

32 

4,5' 



Z2 = l 

Z4 = 1 



Sffi 



Jh 

16(5 8u(l + u)e' 

32 

2(5' 



(23) 



from which and the relations ( 19 1 we can deduce the renormalization constants for 



the fields and parameters of the model 
3.92 , 3gi 



Zr- 
Zd 



1 + 

: 1 



Zg^ — 1 



166 8u(l + u)e' 

3gi g^ 

8u{l + u)e 16(5' 

3.g2 3ffi 

4(5 4M(l + u)e' 



Z. 



z. 



z. 



■4>'i 



92 
16(5' 



Zu = l 



3.92 



32 



3.91 



8u{l + u)e' 
3.91 



16(5 8M(l + u)e 



(24) 



5 Fixed points 

The coefficient functions of the RG operator 



^RG = M 






d_ 



t'l^.+^P^l^^loD 



d_ 
'dg. 



d_ 
do' 



(25) 



where the bare parameters are denoted with the subscript "0", are defined as 



IF = M- 



dlnZp 



d„. 



Pi= 1^ 



(26) 



with the charges gi = {(7i,g2,w}. From this definition and from relations (17 1 
follows the explicit form for the beta functions 

(3g,=gii-2€-r] + S-fD), I3u = u[-rj + -fo], ^g, = g2(-2(5 - 7sJ (27) 
and for the relevant anomalous dimensions 'jjj and 7^2 
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"2u(l + u) 



332 
2 ■ 



(28) 



The scaling regimes are associated with the fixed points of the corresponding RG 
functions. The fixed points g* are defined as such points g* = {g^, m*, (72) for which 
all /3g functions vanish 

/?si (9*1 , w* , ff2 ) = Pu (31 , M* , 32 ) = P92 (9*1 , u* , 32 ) = 0. (29) 

The type of the fixed point is determined by the eigenvalues of the matrix fi = 



{rtik = dPi/dgk}, where /3i is the full set of /S functions (27) and gk is the full set 

of charges {gi,u, 172} • The IR asymptotic behavior is governed by the IR stable 

fixed points, for which all eigenvalues of il matrix are positive. 

It is easy to see that the functions f3g-^ and /3„ satisfy relation Pg^/gi — 3/3„/u = 

2(?7 — e). This means that they cannot be equal zero simultaneously for the finite 

values of the charges gi and u. The only exception is the instance e = rj, which 

should be studied separately. For general case e 7^ 77 we have to set either u = 

or u = 00 and rescale g in such a way, that jd remains finite |20| . 

In what follows we present the results for fixed points, anomalous dimensions and 

eigenvalues of the Jl matrix to the first order of perturbation theory. However, 



we would like to stress, that the form of /3 functions (27 1 allows to calculate the 
anomalous dimensions jjj and 7^ exactly (without any second-order correction). 
In |28| independence of the renormalization constant Zu on the exponents 77 at 
the two-loop approximation has been conjectured. It implies that we may use the 
choice 7] — 0, which we have applied in our calculations of the renormalization 
constants Zd and Zg^ . 

Let us consider the "rapid-change mode" (u -^ cx)). It is convenient to introduce 
new variables w — l/u,g'i — gi/v? and the corresponding j3 functions obtain the 
form 

/^si =3i[7?-2e + 7i3], /3^ = ^[77 - 7d], /^g^ = 32(-25 - 7gJ. (30) 

The "rapid-change model" corresponds to the fixed point with w* = 0. In this case 
four stable IR fixed points are realized: 

FPIA: w*=Q, g[*=0, 9*2 =Q (31) 

Id = 0, 7s2 = 0- 



(32) 



stable for 


ry > 2e, (5 < 0, 77 > 


FP IB: 


w*^Q, 5r=0, 9*2 = - 




715 =g, ig^=-25. 


stable for 


677 -f (5 > 12e, (5 > 0, 5>Qr 


FP2A: 


»-o. A'^''\-'\ 




-fD = '2e- 77, 7g, = 277 - 4e 


stable for 


2e > 77 > e, 2e> rj + S 



9*2 = (33) 
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PP.B: „,.^0. ,;• . fflHi_|!l^, ,;.ffi±l^ (34, 

15 5 

7r, = 2e - -q, jg^ = -2A 

stable for S > 0, S + ij > 2e, 12e > 6t] + S 

For the analysis of the regime w — > (quenched velocity field) we introduce the 
new variable gi" = gi/u. Hence the corresponding (3 functions have the form 

/^si" =5i"[-2e + 27£,], /3„ = ^[-77 + 7^], I3g^ = g2{-2S - •jg^) (35) 

Also in this case there are four possible IR stable fixed points: 

FP 3A: u* = 0, g'l* =0, ^2 = (36) 

ID = 0, 7g2 = 0. 



stable for 


e<0, 


(5<0, 


7/ <0 


FP3B: 


u* =0, 


11* 
.91 = 


0. ».^f 




6 
7D=g 


> 7s2 = 


= -2A 


stable for 


e< 0, 


(5>0, 


r/<0 


FP 4A: 


w* =0, 


.9i == 


^|, »^^o 




Id = e, 


7g2 = 


-2e 


stable for 


e> 0, 


6X5, 


e>r/ 


FP 4B: 


w* =0, 


3i = 


4(6e - (5) 
15 




Id = e, 


7g2 = 


-2A 


stable for 


0< e< 


(5<6e, 


e > 77 



(37) 



(38) 



52 = s (39) 



In the special case e = 77 the functions f3g-^ and /3,j become proportional and 
this leads to the degeneration of fixed point. Instead of just plain fixed point, we 
observe a whole line of fixed points in the {gi,u) plane. 

FP5A: f[ = ^, g;=0 (40) 

u*[l + u*) 6 

715 = e = 77, 7g, = -2e 
stable for r] = e > S, e>0 

FP^B: ^v^l^^i^, *^«^ (41) 

u*[l + u*) 15 5 

7-D = e, 732 = -2(5 
stable for e ^ rj < S, e > 0. 

The "real problem" corresponding to the Kolmogorov scaling is obtained for 
the value rj = e = 4/3, which leads to the famous "five-thirds law" [27^ for the 
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spatial velocity statistics. By direct observation we see that in this case critical 
behaviour is described by the fixed point 5A for logarithmic((5 = 0,d = 4), three 
dimensional ((5 = 1/2, d — 3) and also two dimensional(5 — l,d — 2) case. Fixed 
point 5A is characterised by vanishing of the charge g2. Therefore we conclude 
that in the turbulent field cubic interactions in the vicinity of critical point (|3| are 
in fact negligible. 



6 Conclusions 

This paper is devoted to the study of directed percolation problem influenced by 
the external advecting velocity field. In order to use the technique of the pertur- 
bative renormalization a field-theoretic model is constructed. All the calculations 
were performed to the first order of the perturbation theory. The IR stable fixed 
points, dimensions and corresponding regions of stability of fixed points are cal- 
culated. The technically relatively simple model of velocity fluctuations used here 
is a convenient starting point for more realistic high-loop calculations. 
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